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Abstract. We review interaction models for GRB afterglows, with an eye on constrain- 
ing the nature of their progenitors and the geometry of explosion. Evidence is presented 
for two types of progenitors. The radio afterglow of GRB 980425/SN 1998bw and the 
multi-frequency observations of GRB 970508 can be fitted by a blast wave expanding 
into a wind-type medium, pointing to a massive star progenitor. The broadband af- 
terglow data of GRBs 990123 and 990510 are better modeled by a jet expanding into 
a constant-density medium, implying a compact star merger origin. Among other well 
observed GRBs, the jet model appears to be more widely applicable than the wind 
model, although some cases are ambiguous. The model fits often require a deviation 
of the energy distribution of the radiating electrons from the commonly assumed sin- 
gle power-law form, particularly for sources with a rapid decline and/or pronounced 
steepening in the optical light curves. Transition to non-relativistic evolution has been 
suggested as an alternative explanation for the light curve steepening, although to pro- 
duce steepening on the order of days or less would require very high ambient densities, 
which are generally difficult to reconcile with radio observations. Major open issues 
include the hydrodynamics of jet-ambient medium interaction, a self-consistent deter- 
mination of the electron energy distribution, and the effects of pair production on the 
early afterglows which are expected to be particularly large for the wind interaction 
model. 



1 Introduction 

The simplest model of the afterglows of gamma-ray bursts (GRBs) involves a 
spherical relativistic blast wave expanding into a constant-density, presumably 
interstellar, medium [ ^8[|7l[ |. The afterglows are emitted by nonthermal electrons 
accelerated at the shock front to an energy distribution usually assumed to be 
a power law above some cutoff determined by the shock velocity. The predicted 
power-law decay of the af terg low emission with time was subsequently observed 



at X-ray |gO[ and optical |113| waveleng ths for G RB 970228, giving basic confir 



mation to this now "standard" picture |116;125 



To date, nearly 50 GRB afterglows have been observed at more than one 
frequency (for reviews of afterglow observations see Chapters by Piro, Plan & 
Fruchter, and Weiler). A dozen or so of these are well observed at multiple fre- 
quencies to allow for detailed modeling. Some features of these afterglows turn 
out to be difficult to accommodate in the standard model. The most noticeable 
is the relatively steep decline of the optical light curves, sometimes preceded by 
flatter evolution. The steep decline is usually attributed to a collimated or jet- 
like, rather than spherical, initial energy injection ^ ,105|. Alternatively, it could 



2 Li and Chevalier 



be due to the transition of blast wave evolution to the non-relativistic regime 
[ p6| , a non-standard electron energy distribution [ |66||7^ , and/or expansion into 
a wind-type ambient medium p^ . The question of a wind versus a constant- 
density surrounding medium is a crucial one for the progenitors of GRBs, since 
massive stars, one of the leading candidates for GRB progenitors [127 7^, should 
be surrounded by a wind. In contrast, GRBs resulted from compact star merg- 
ers, the other leading candidate, are expected to be surrounded by the general 
interstellar medium (see Chapter by Waxman). 

The plan of the Chapter is as follows. In § ||, we discuss the likely environ- 
ments of GRBs relevant to afterglow evolution. In § |^, we outline the basics of 
various ambient interaction models for GRB afterglows. These models are then 
applied to individual GRBs with relatively well-observed afterglows (§ ^). We 
find that some of the sources are probably wind-interactors, while others are 
better modeled as interacting with a constant-density medium. This and other 
results are discussed in § |[ 



2 The Ambient Medium 

Models for the GRB afterglows indicate that the emission comes from a region 
~ 10^^ — 10^^ cm from the source of the explosion. The nature of the material 
in this region depends on the GRB progenitors, which are presently not known. 
There are two main possibilities. In one model, the progenitors are massive stars 
at the end of their lives [ 
primarily with the mass loss from the progenitor leading up to the explosion. 
The other class of object is the merger of two compact objects for which there 
are a number of possibilities |3^. In this case, there may be some debris from the 
merger, but the interaction is expected to be primarily with the ISM (interstellar 
medium). We consider these possibilities in turn. 

If the progenitors are massive stars, there is an analogy to the explosions 
of core collapse supernovae, for which there is abundant evidence that they are 
interacting with the winds from the progenitor stars (see Chapter by Chevalier 
& Fransson). In general, the interaction appears to be with the free wind from 
the progenitor star, with density pw — M / {'inr'^Vw), where M is the mass loss 
rate from the star and Vw is the wind velocity. In most of the supernova cases, 
the radial range that is observed is out to a few 10^^ cm, so that the mass 
loss characteristics have not substantially changed during the time that mass is 
supplied to the wind and pw oc r~^. The density in the wind depends on the type 
of progenitor. Red supergiant stars, which are thought to be the progenitors of 
most Type II supernovae, have slow dense winds. Wolf-Rayet stars, which are 
thought to be the progenitors of Type lb and Ic supernovae, have faster, lower 
density winds. SN 1987A is a special case in which the star was a red supergiant, 
but became a blue, B3 supergiant about 10^ years before the explosion (see 
Chapter by McCray). The result was a complex circumstellar medium, including 
a dense ring at a radius of 6 x 10^^ cm. 



127:73]. In this case, the interaction is expected to be 
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If GRBs do have massive star progenitors, there are a number of arguments 
suggesting that Wolf-Rayet stars are the most likely progenitors: (1) SN 1998bw, 
the best case of a SN - GRB association (GRB 980425), was of Type Ic, with 
a probable Wolf-Rayet progenitor. (2) The high energy of GRBs suggest that a 
moderately massive black hole is involved, which, in turn, requires a massive, 
^20 — 25 Mq, progenitor These stars are likely to be Wolf-Rayet stars at 
the end of their lives |Q. (3) The relativistic flow from a central object may 
be able to penetrate a relatively compact Wolf-Rayet star, but probably cannot 
penetrate an extended red supergiant star ||6|]. (4) If a rapidly rotating black 
hole is required for the GRB explosion, a merger at the center of a massive star 
is one way to obtain a rapidly rotating core. The merger process can give rise 
to a Wolf-Rayet star. None of these arguments is definitive, but they do point 
toward the most plausible progenitors. 

The circumstellar medium created by a 35 Mq initial mass star that passes 
through a Wolf-Rayet phase and a prior red supergiant phase has been simulated 
by Garcia-Segura et al. Q. They find that during the Wolf-Rayet phase, a 
broken shell is created at a radius of a few pc or more. This type of shell has 
been observed around a number of Wolf-Rayet stars. Inside of the shell is a 
region of shocked wind, with an approximately constant density The size 
of the shocked region is determined by what is required to decelerate the wind 
and is typically less than half of the shell radius. In general, the wind ram 
pressure p — PwV^ determines the size of the wind bubble that is created. The 
pressure is p/k = 3.6 x 10^ M-^v^r^^ cm"^ K, where M_5 = M/10"^ Mq yr"\ 
Vs, = Vw/W^ cm s""'^, and riy ~ r/lO^^ cm. For comparison, the interstellar 
pressure in the solar neighborhood is ^ 3000 cm^"^ K. The region out to the 
wind termination shock can be described by a density pu, = 5 x 10^^A*r~^ g 
cm~^, where = IM-^v^^. Considering the radial range relevant to GRB 
afterglows, we expect this to apply to most cases, unless the wind is especially 
weak or the surrounding pressure is extraordinarily high. Ramirez-Ruiz et al. 

have followed the circumstellar mass loss from Wolf-Rayet stars and found 
that the free wind does not extend to large radii in their treatment. 

For bursts that interact directly with the ISM, the size scale of the afterglows 
suggests that the medium has a roughly constant density, although the ISM is 
known to have considerable density inhomogeneity on a range of scales. The 
range of interstellar densities in a galaxy like our own goes from ~ 10"'^ cm~^ 
in the hot ISM to >10^ cm~'^ in compact cores in molecular clouds. The high 
densities are found in a very small volume fraction so that a typical event is 
unlikely to occur in such a region. The situation may change in a starburst 
region, where a significant fraction of the gas may have a density <;10'^ cm~'^ 
[ |107| . The radio-bright, compact supernova remnants in starburst regions may 
be interacting with this dense interstellar component jl^ . The finding that GRB 
010222 may have occurred within a very active starburst indicates that the 
interstellar environments of GRBs may have unusual properties. 

Although the interstellar and wind models are the two main types of environ- 
ments considered for afterglows, there is a different scenario involving a massive 
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star, motivated by the observation of possible Fe lines in the X-ray spectra of 
afterglows 128 8^,|^. The observations appear to require a substantial mass 
of Fe at r 10^ cm, so one suggestion is that a supernova occurs before a 
GRB and that the ejecta have time to expand to this radius |114 . This requires 
that the GRB be delayed by months. The supernova would expand into the 
progenitor wind, creating a complex circumburst region in the inner part of the 
wind. 



3 Ambient Interaction Models 
3.1 Standard ISM-Interaction Model 

We start with the simplest model of GRB afterglows involving synchrotron emis- 
sion from a spherical relativistic blast wave propagating into a constant-density 
interstellar medium. The model was worked out in some detail by Meszaros & 
Rees |7l| in advance of the afterglow detection. It has since been elaborated 



upon by others (e.g., [116 117,124|), and in particular 106 . Here, we will follow 
the formalism of Sari et al. |l06| , which has been widely used to interpret the 
spectra and light curves of GRB afterglows. This standard model serves as a 
benchmark against which other models will be compared. 

The standard model assumes that the synchrotron-emitting electrons are 
accelerated at the shock front to a power-law distribution of Lorentz factor 
iV(7c) djc oc (with p > 2) above some minimum cutoff 7ni, which 

is determined by the shock velocity. This distribution is further modified by 
synchrotron cooling in the downstream flow. Denote by 7c the Lorentz factor of 
an electron which cools in a blast wave expansion time. In the "fast cooling" case 
with 7c < 7m, the electron distribution is given by 7,7^ between 7c and 7^ and 
7(7''^^ above 7m- This electron distribution produces a synchrotron spectrum of 
the flux F,y oc between fc and fm, the characteristic frequencies of the 

photons produced by electrons with 7c and 7m respectively, F,y oc v~'pI'^ above 
fc'm, and Fi, oc v^l'^ below v^- The spectrum is modified by synchrotron self- 
absorption, which typically occurs at relatively low frequencies. Assuming that 
the self-absorption frequency is less than the cooling frequency and that the 
absorbing electrons that have cooled for different lengths of time are well mixed 
spatially |Q, we have Fy oc below v^. In the opposite, "slow cooling" case 
with 7c > 7m, the electron distribution steepens to 7^^"^ above 7c. Assuming 
that the self-absorption frequency is less than the characteristic frequency I'm, 
we have the following spectrum: F^, oc for v < z^a, F^, oc z/^/'^ for i/g, < v < fm, 
Fi, oc for z/m < I' < t'c, and F^, oc v~'pI'^ for v > v^- These broken power- 

law expressions are often used to interpret the instantaneous spectra of GRB 
afterglows. They tend to agree with the more accurately determined spectra to 
within a factor of a few (e.g., [^ ). 

To obtain light curves, one must first determine the time evolution of the 
characteristic frequencies (t'a, and v^) and the peak flux i^i/,max- The evolu- 
tion depends on the dynamics of the blast wave. For the spherical blast wave 
propagating into a constant-density medium envisioned in the standard model, 
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the dynamics are described by the self-similar solutions of Blandford & McKee 
[|j : the blast wave Lorentz factor 7 varies with the distance R from the explosion 
center as 7 cx R~^^^ in the adiabatic regime and 7 oc in the fully radiative 
regime. The Lorentz factor and distance are related to the detector-time t as 
t oc Rjy^ ■ Making the standard assumption that a constant fraction (cb) of 
the blast wave energy goes into the electrons (magnetic fields) and making use 
of the jump conditions for relativistic shocks, one obtains the following scalings: 
z^a oc oc OC t~^l'^, and F^,^^-^ oc t° for the slow cooling case, 

which is always in the adiabatic regime. In the fast cooling case, if the electron 
energy fraction Cc is close to unity, then the blast wave is fully radiative, and 
the scalings become: i^a cx t^^l^ , oc t"^/'', oc t~^'^l'^ , and F^^niax oc t~^^'^ . 
If on the other hand eo ^ 1, the evolution remains adiabatic, and the scalings 
are the same as in the slow cooling case, except for the self-absorption frequency 
which now decreases with time as fa oc t~^/'^ . These scalings are combined with 
the instantaneous spectra listed above to obtain light curves. 

The light curve at a given observing frequency v can be described by a broken 
power-law with breaks at various characteristic times. These include the times 
ta, im and tc when the characteristic frequencies Vg., and pass, respectively, 
the observing frequency v. Ignoring self-absorption. Sari et al. [106| obtained two 
types of light curves in two frequency regimes, separated by a critical frequency 
Vq at which the characteristic times t-^ and become equal, i.e., <m(t'o) = 
tc{vQ) = to. The critical time to divides the early, fast cooling part of the light 
curve from the later, slow cooling part. In the high frequency regime where 
u > uq, the ordering of the characteristic times is ic < < ^Oi which yields a 
light curve of F,, oc t^/^ for t < t^, oc t-^/'^ for < t < t„„ and F^ oc t(2-3p)/4 
for t > tm 'vn. the case of adiabatic evolution. In the opposite case of fully radiative 
evolution, which is possible only before the critical time when the electrons are 
fast cooling, one has instead oc t^^/^ for t < tc, F^, oc i"'*/^ for tc < t < im, 
and Fi, oc ^(^-ep)/? Jqj, < ^ < jn the low frequency regime where 1/ < vq, 
the ordering of the characteristic times becomes t^ < t^ < tc- The light curve 
is then oc t^/^ for t < to, F^ oc t^/^ for to < t < t^, F^ oc i3(i-p)/4 for 

tra <t < tc, and F^ oc t(2-3p)/4 fQj. ^ 

> tc in the adiabatic case. In the fully 
radiative case, one has F^^ oc t^^/-^ for t < to instead. 

Synchrotron self-absorption modifies the light curves. We will concentrate 
on the case of adiabatic evolution here and below; it is more likely than the 
fully radiative case during most of the afterglow phase of a GRB evolution (e.g., 
[[71]). The self-absorption is characterized by the absorption frequency fa, which 
decreases with time as fa oc t~^/^ during the fast cooling period t < to. It has 
a constant value, denoted by fa,0: between tg < t < ^am, where ^am is another 
critical time when the characteristic frequencies fa and fm become equal. In the 
frequency regime f > fa,0: both the high and low frequency light curves discussed 
above have an additional power-law segment F^ t before the characteristic 
time tg_, which is smaller than any other characteristic time. Collecting all time 
and frequency dependences, we finally have Fi, oc tv^ {t < ia), F^, oc t^^^v^^^ 

{t^<t< tc), F^ oc (tc<t< tm), and F^ oc t(2-3p)/4j,-p/2 > 
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for the high frequency hght curve, and F^, oc tv^ {t < ta), F^, oc t^/^v^^-^ (t^ < 

t < to), F, oc (^to <t < t„), F, oc t3(l-p)/4^-(p-l)/2 (^^^ < ^ < 

and oc t(2~3p)/4^-p/2 ^-^ ^^-j j^j, frequency light curve. The time and 

frequency dependences of the flux are sometimes parameterized using F^ ozf^v^ . 
The relations between the decay index a and spectral index (3, a = 3(1 — p)/4 = 
3/3/2 (t,„ <t <tc or < V < ly^) and a = (2 - 3p)/4 = (1 + 3/3)/2 (i > or 
ly > i^c > i^m) in the most relevant, slow cooling case, are often used to interpret 
optical and X-ray afterglow observations at relatively late times (of order one day 
or more). In the frequency regime v < v^.o, more relevant to radio afterglows, 
one has the ordering to < tm < < tc, which yields F^ oc iz^^ {t < to), 

F, oc tl/2j,2 (^to<t< t^), F, oc ^5/41.5/2 (<,„ < t < t^), F, OC t3(l-p)/4j,-(p-l)/2 

{t^ <t < tc), and F^ oc t(2-3p)/4j,-p/2 > reiterate that the scahngs 

involving self-absorption in the fast cooling case with t < to are derived assuming 
that the energy distribution of the cooling electrons is spatially homogeneous in 
the emission region behind the shock front. The opposite situation where layers 
of cooling electrons at different distances from the shock front remain unmixed 
has been considered in detail by Granot et al. [|8| and Granot & Sari pof . 

3.2 Wind-Interaction Model 

The first modification of the standard model we consider is the density distribu- 
tion of the ambient medium. Meszaros, Rees, & Wijers |Q studied the general 
case of a power- law ambient density distribution n oc r~*, with an arbitrary 
power-index s. Chevalier & Li |T7|[r^ examined the specific case of s = 2, cor- 
responding to a constant-mass loss rate, constant-velocity, circumstellar wind, 
possibly of a Wolf-Rayet origin. Some features of the wind-interaction model are 
also described in Dai & Lu |^ and Panaitescu, Meszaros & Rees |84j. Here, we 
follow the formalism of Chevalier & Li [ p^ . 

The wind-interaction model has the same instantaneous afterglow spectra 
as in the standard ISM-interaction model, but different light curves. The differ- 
ences in light curve come from the blast wave dynamics. For a relativistic blast 
wave propagating in an s = 2 medium, Blandford & McKee showed that 
its Lorentz factor 7 oc R~^^^ in the adiabatic regime, and 7 oc in the fully 
radiative regime. Concentrating on the adiabatic evolution as before and making 
the standard assumption about the fractions of the blast wave energy going into 
the electrons and magnetic fields, one finds that the characteristic frequencies 
scale with time as oc t^^^^, v^^ t^"^/^, and i>c oc t^^^, and the peak flux 
-fl/.max oc t~^/2 slow cooling case. In the fast cooling case, the scalings 

are the same except for the self-absorption frequency, which now decreases with 
time more rapidly as i^a oc t~^/^. 

The above scalings of the characteristic frequencies define four critical times, 
with four corresponding critical frequencies which divide the light curves into five 
distinct frequency regimes. For the typical parameters adopted by Chevalier & 
Li [|l8| , the light curve that is most relevant to the optical and X-ray afterglows 
has the ordering ta < < to < tc, which yields F^, oc tJl'^v^l'^ for t < t^, 

F^ oc t-l/4^-1/2 foj, ta < t < tm, F^ OC t(2-3p)/4^-p/2 foj. ^^^^ < t < to, F^ <X 
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^(2-3p)/4j^-p/2 fQj. ^ t < tc, and OC i(l-3p)/4^-(p-l)/2 foj. ^ > 

that as the cooling frequency Vc moves across the observing frequency v from 
below, the Hght curve steepens from a = (2 — 3p)/4 to (1 — 3p)/4 by a modest 
amount Aa = 1/4. The same amount of steepening, from a = 3(1 — p)/4 to 
(2 — 3p)/4, also occurs in the standard ISM-interaction model, as I'c moves 
across v from above. A signature of the wind-interaction model is the relatively 
fast decline of the light curve in the non-cooling spectral region ly^ < ly < I'c, 
where a = (1 — 3p)/4 compared with 3(1 — p)/4 for the ISM case. If the cooling 
frequency falls between optical and X-ray wavelengths, the wind model would 
predict an X-ray light curve declining less steeply than an optical light curve, 
whereas the opposite would be true for the ISM model. 

At the much lower, radio frequency, the characteristic times typically have 
the following ordering: < < < ^m, which leads to a light curve with 
Fi, OC t'l^v^l'^ for t < tc, Fy OC t^v'^ for < t < to, F^ oc tv'^ for to < t < tg,, 
F^ OC t°iy^/^ for < t < tn,, and F^ oc i(i-3p)/4j,-(p-i)/2 j^^. ^ > ^j^g 

slow cooling case with t > to most relevant to radio observations, the fairly 
steep rise of radio flux F^, (x t at the self-absorbed frequencies could in principle 
distinguish the wind model from the standard ISM model, where the rise is 
slower {Fi, oc t^^^). In practice, the difference is masked to a large extent by 
interstellar scintillation. 



3.3 Transition to Non-relativistic Evolution 

After a GRB blast wave sweeps up a mass equivalent to the rest mass of the 
explosion, its evolution becomes non-relativistic. In a medium of constant density 
n, this happens at a r adius of order R^r- = 1.2 x IO^^eI^^ 'n cm where £'52 
is the explosion energy in units of 10^^ ergs and n the ambient number density 
in units of cm""^. The transition occurs on a time scale t^r = Rm/c = 1.2(1 -|- 
z)i?52^ri~^/^ yr (where z is the cosmological redshift), which is longer than the 
typical duration of GRB afterglow observations, unless the explosion energy is 
much lower than lO^'^ ergs |Q and/or the density much higher than 1 cm~^ p^ . 
Provided that the minimum electron Lorentz factor 7,„ > 1, the temporal and 
frequency dependences of the flux in the non-relativistic regime can be obtained 



in a way similar to that in the relativistic regime [125|. Here, wc follow the 
detailed treatment of Frail et al. |3^. The electrons arc typically slow cooling, 
with t'm oc t~^ , fc oc and the peak flux at v^^ given by i^z^.max k t'^^^. 

The self- absorption frequency could either be greater or smaller than i^ni- There 
are two relevant orderings of characteristic frequencies: I'a < i^m < i^c and v^^ < 
Va. < i^c- In the former case, i^a i^/^, one has oc v^t'"^!^ for v < i^a, 

F^ oc 1/1/3^8/5 for < y^^ Fy oc jy-(P-l)/2i3(7-5p)/10 foj. j,^ < < jy^^ and 

Fy (X V 

-p/2^(4-3p)/2 foj. ^ ^Yie latter case, l^a oc t^'^ 3p)/(p+4)^ g^j^j ^ 

1.2^13/5 for ^ < oc ^.5/2^11/10 for u^<v < l/a, OC j.-(p-l)/2t3(7-5p)/10 

for Vg. < V < Vc, and F^ oc iy^p/2i(4-3p)/2 foj. ^ Note that the hght curve 
decreases with time faster in the non-relativistic regime (where a = 3[7 — 5p]/10 
and [4 — 3p]/2 before and after the cooling break t^ ot a — —1.65 and —1.75 for 
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p — 2.5) than in the relativistic regime (where the corresponding a = 3[1 — p]/4 
and [2 - 3p]/4 or a = -1.125 and -1.375). 

In the wind interaction model with an ambient density p — Ar^'^, the transi- 
tion to non-relativistic evolution occurs at a distance i?nr = 1-8 x 10^*£^52/^* cm, 
where is the coefficient A in units of 5 x 10^^ g cm~^ [Q, corresponding to a 
time tni = 1.9{l + z)E52/At, yr, again longer than most afterglow observations for 
typical parameters. In the non-relativistic regime, one finds oc t^^/^, i^c oc i, 
and the peak flux Fj^,max oc t~^^^ at Vrn- The relevant orderings are i^a. < i^m < Vc 
and v-ax < Va < Vc- In the former case, Va, oc and oc v^t ^ for V < I^a, 

oc 1/1/3^4/9 foj. j,^ < J, < jy^^ oc i.-(P-l)/2i(5-7p)/6 foj. < < j,^^ and 

i^i, oc j,-p/2t(8-7p)/6 foj, ^ > j^^^ ^Yie latter case, oc t(4-7p)/[3(p-K4)] ^ ^^^^^ 

i^^ oc vH-^ for 1/ < t/ni, F^ oc iy5/2tl/6 for < < ly^^ oc ,y-(p-l)/2i(5-7p)/6 

for < < Vc, and oc p-p/'^t'^^-'^p)!^ for > i^c- As in ISM interaction, 
the light curve decreases with time faster in the non-relativistic regime (where 
Of = [8 — 7p]/6 and [5 — 7p]/6 before and after the cooling break tc, or a = —1.58 
and —2.08 for p = 2.5) than in the relativistic regime (where the corresponding 
a = [2 - 3p]/4 and [1 - 3p]/4 or a = -1.375 and -1.625). 



3.4 Jet Model 

Starting with Rhoads ||9^], jet models of afterglows have been widely discussed in 
connection with the steepening of light curves. We will first outline the asymp- 



totic analysis of Sari, Piran & Halpern |105 and then comment on possible 
complications. Recent reviews of the subject include and 

Let 9o be the initial angular width of the jet. When the Lorcntz factor 7 
drops below Oq^, the jet starts to spread sideways, changing its dynamics. The 
spreading occurs around a time tjct — 3[(l-l-2:)/2](£'52/n)^/"^(6'o/0.2)*/'^ days in a 
constant-density medium and ijot — 2[(l-f z)/2](_E52/A,)(0o/O.2)^ days in a wind 
[ p^ . Before ijet, the standard spherical results apply. After tj^t, the jet Lorentz 
factor decreases with distance exponentially |Q, yielding 7 oc z/m oc 

and Vc oc t", and the peak fiux i^,y.max t^^ at v^. In the case < t'm < ^'c, 
one finds that i/s. oc t~^/^, and Fi, oc ly^t^ for ly < v^, F^ oc v^/^t^^/^ for 
Va<v< v-m, Fi, oc i/-(P-i)/2t-P for <v < Vc, and F^ oc v^P/'^t~P for v > v^- 
In the case t'm < < Vc-, we have instead oc t^2(p+i)/(p+4) ^ ^-^^ ^ ^2^0 

for V < Vra, Fy oc V^/'^t for V^a < V < Va,T Fy (X :y-(P-l)/2i-P for I/a < < Vc, 

and oc v^P^'^t^P for > z^c. These scalings apply to both ISM and wind 
interaction models. 

The steepening of light curves to an asymptotic scaling of Fy oc is a 
signature of the jet models. How sharply the change from one temporal slope to 
another actually happens remains controversial. Meszaros & Rees |7^ pointed 
out that seeing the edge of a jet when its Lorentz factor drops below 9^^ would 
also steepen the light curve, by a factor of Aa = 3/4 in a constant-density 
medium and Aa = 1/2 in a wind. Semi-analytic calculations taking into account 
of both the sideways spreading and edge effect find that the transition from one 
slope to another tends to be continuous, spanning one decade or more in the 



observer's time |7|,|l2^,||, especially in a wind-type ambient medium. These 
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calculations adopted a simplified set of ID equations for the jet dynamics, which 
need to be checked against numerical simulations. Preliminary 2D calculations 



of Granot et al. |50 show that the blast wave has an egg- like shape, differing 
considerably from that predicted from the analytic or semi-analytic models. A 
steepening of light curves does occur around the time j (x 9^^ , as expected from 
simple arguments. Whether the sharpness and amount of steepening match the 
predictions of the semi-analytic models remain to be seen. 



3.5 Other Modifications 

Many additional effects have been considered in the afterglow literature on top 
of the basic interaction models outlined above. We will limit ourselves to three of 
the more commonly discussed ones: inverse Compton scattering, pair production, 
and a non-standard electron energy distribution. 

Inverse Comptonization of synchrotron photons by relativistic electrons could 
not only produce a high energy component of afterglow emission that is poten- 
tially observable b ut al so dominate the co oling of electrons and affect the blast 



wave dynamics ||83,119, 19y79[| . Sari & Esin |102] examined these effects in detail, 
and concluded that the inverse Compton spectra broadly resemble the primary 
synchrotron spectra in shape, although significant differences do exist, especially 
at the high frequency end where a broken power-law description is no longer ad- 
equate. They showed that as long as the fraction of the blast wave energy in 
electrons e,, exceeds that in magnetic fields ee, the inverse Compton emission 
dominates the synchrotron emission in cooling the electrons in the fast cooling 
regime. Depending on the ratio Sc/^b, the domination can extend well into the 
slow cooling regime, changing the value of the cooling frequency i^c and thus the 
light curves of afterglow emission. 

Thompson & Madau [ |112| considered pair formation as the gamma-ray pho- 
tons from the GRB proper propagate ahead of the external shock front and 
interact with the seed photons back-scattered by the ambient medium. A sim- 
plified discussion of the process is presented in Meszaros, Ramirez-Ruiz & Rees 
[ [70| . Beloborodov Q went one step further and solved for the dynamics of the 
pair-loaded medium. These studies concluded that pair production can dramat- 
ically increase the radiative efficiency of the blast wave and potentially broaden 
the original pulses of gamma-rays if the ambient density is high enough. The 
effects on afterglows are expected to be large at early times, especially for the 
wind-interaction model, although details are yet to be worked out. 

Most afterglow models assume that the radiating electrons have a power- 
law energy distribution at injection above some cutoff. The power-law indexes 
inferred from afterglow observations under this assumption span a wide range, 
from p ~ 2 (or below) to 3 (e.g., |]l8| , ^0|j8ll| ) . The lack of a universal value for 
p calls into question the assumption that the shock front accelerates electrons 
to a power-law that is constant with energy and with time. Indeed, in the Crab 
Nebula, arguably the best studied astrophysical synchrotron source, a break 
in the injection spectrum is required ||l|. This motivated Li & Chevalier |]66| 
to consider a broken power-law distribution of electron energy that steepens at 
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high energies in the context of the spherical wind-interaction model. Similar non- 
standard energy distributions have been adopted in other models [24,78 12l| . A 
major attraction of the invoked steepening of electron energy distribution is that 
it can lead to a steepening of afterglow light curves. The light curve steepening 
is chromatic, and should be distinguishable from the achromatic steepening due 
to jet effects if wide wavelength coverage is available. 



4 Application to Individual Sources 

We now apply the interaction models to the dozen or so GRBs whose afterglows 
are relatively well observed. It turns out that for most of the sources the models 
are not unique. We shall start the discussion with the cases that in our view are 
less controversial. 



4.1 Probable Collimated ISM Interactors: GRBs 990123 and 990510 

The first two GRBs that show breaks in their optical afterglows are GRB 990123 
[ IsT] and GRB 990510 ||5^. The breaks in both sources are achromatic and are 
interpreted as due to jets. In the case of GRB 990123, the break occurs at a 
time ti, = 1.68 ± 0.19 days, depending somewhat on the fitting function adopted 
Before the break, the temporal decay index at optical wavelengths is a = 
-1.12 ± 0.08. The optical spectral index is found to be /? = -0.750 ± 0.068. 
These two indexes can be fitted by an ISM- interaction model with p = 2.5 in the 
adiabatic regime (where a = ~3[p — l]/4 and /? = — b — l]/2 for i^m < < i^c). 
The decay index a = — 1.44±0.07 at X-ray wavelengths is also consistent with the 
model, provided that the X-rays are in the cooling regime (where a = [2 — 3p]/4 
for V > i^c)- The observed steeper decline in X-rays than in optical is expected for 
ISM interaction but not for wind interaction ||l^. After the break, the optical 
light curves steepen quickly to a = —1.69 ± 0.06 while the spectral index (3 
remains approximately the same The observed amount of steepening, Aa = 
0.57 ± 0.10, agrees marginally with that expected of a jet with a fixed opening 
angle (where Aa = 3/4 for ISM interaction) . It is less than half of the asymptotic 
value Aa = (p + 3)/4 = 1.38 expected from the sideways expansion of a jet. The 
relatively small amount of steepening was taken as evidence for seeing the edge 
of a non-spreading jet |j7^ , ^l|j5^ ] , although it does not rule out jet-spreading as 
the cause of steepening, since it takes time for the light curves to reach the 
asymptotic slope. Indeed, Panaitescu & Kumar obtained a reasonable fit to 
the multi-frequency data of GRB 990123 using a semi-analytic jet model taking 
into account of lateral expansion. 

The breaks in the optical light curves of GRB 990510 appear to be smoother 
than those in GRB 990123. Fitting a continuous function, Harrison et al. |Q 
obtained a break time of ib = 1-20 ± 0.08 days. Before the break, the decay 
index a = —0.82 ± 0.02, which implies p = 2.1 in the adiabatic regime for 
ISM-interaction. The expected spectral index /? — —0.55 is consistent with /3 = 



0.61±0.12 determined by Staneket al. [110| or /3 = -0.531±0.019 by Holland et 
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al. |£5[. The latter authors find some evidence that (i decreases with time which, 
if real, is not explained. The relatively flat decay rate of a > —1 is difficult 
to accommodate in the wind-interaction model; it requires a very flat electron 
energy distribution with p < 2 JlSt . Even then, the predicted spectral index 
would not match that observed llSOl. After the break, the measured a = —2.18 ± 



0.05, which is consistent with the expected asymptotic value a — —p |105|. 
The ISM-interacting jet model is further supported by radio data, which are 
consistent with the expected oc t~^/^ evolution X-ray data are available 
between 0.3 and 2 days, bracketing the break in optical [Q. The lightcurve 
can be fitted either by single power-law with a — — 1.42 ± 0.07 or a broken- 
power law, and is well modeled using an ISM-interacting jet model with the 
cooling frequency Vc between the optical and X-rays pq,p5|. The detection of 



polarization from this event [21 126 is also consistent with the jet interpretation 



[ [45yi02|| , although the observed polarization is small and does not require a jet. 

4.2 Probable Spherical Wind Interactors: GRBs 980425 and 970508 

The probable association of GRB 980425 with the radio supernova SN 1998bw 
at z = 0.0085 is discussed in the Chapters by Galama and Iwamoto & Nomoto. 
Here, we are concerned with the radio emission from this source, which is ex- 
traordinary Kulkarni et al. argued that the shock responsible for the 
radio emission is relativistic, based on the high brightness temperature, on a 
synchrotron self-absorption interpretation of the early evolution, and on scintil- 
lation results. Using the standard synchrotron theory, Li & Chevalier were 
able to deduce from the radio data that the shock is expanding into a medium 
with an approximately density profile at a speed comparable to the speed 
of light. A more detailed model, taking into account of the shock dynamics and 
relativistic effects and assuming the standard power-law electron distribution 
above some cutoff, reproduces the radio lightcurves after about day 10 reason- 
ably well, provided that the shock energy is increased by a factor of ^ 2.5 at 
~ 10^ days in the rest frame of the explosion to explain the rise in the radio 
fluxes observed between days 20-40. Based on the radio spectra at day 12 and 
day 15 and adopting a mono-energetic distribution for the synchrotron-emitting 



electrons, Waxman & Loeb [118| inferred a shock speed of ^ 0.3 c, which is a 
factor of ~ 2 smaller than that inferred by Li & Chevalier around the same 
time. The longer term evolution of the radio source implies a power law particle 
energy distribution, as generally observed in GRB afterglows. Weiler, Panagia & 



Montes |122| were able to fit a parameterized model to the radio data, and came 
to a similar conclusion that the shock is mildly relativistic and is interacting 
with a clumpy progenitor wind. 

Chevalier & Li [|l8| applied the spherical wind-interaction model to GRB 
970508 which, unlike GRB 980425, is at a cosmological distance oi z — 0.835 
[ [75) . Extensive optical data up to hundreds of days are available for this source 
(see Chapter by Plan & Fruchter). They follow a power-law after about day 2 
with a = -1.141 ± 0.014 and (3 = -1.11 ± 0.06 The X-rays have a simi- 
lar decay index of a = —1.1 ± 0.1 and the optical/X-ray spectral index is 
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consistent with (3 = —1.1 Q]. Based on the radio to X-ray spectrum on day 
12.1, Galama at ah |^ deduced that the coohng frequency is just below optical 
frequencies. At the optical and X-ray frequencies above the cooling frequency 
V > Vc, the observed value of a implies p = 2.2, which yields a spectral index 
/? consistent with that observed for both the ISM and wind interaction models. 
Radio observations are crucial for distinguishing between these two models. In 
Fig. 1(a), we show the predictions of a wind interaction model at three wave- 
lengths. They fit the radio data reasonably well, especially after about day 100 
when the interstellar scintillation dies down In Fig. 1(b), we present an an- 
alytic fit to the R-band data, which is used in conjunction with the radio fit to 
obtain the best model parameters: eo = 0.2, ee = 0.1, E52 = 0.3, and A* = 0.3. 
The inferred mass loss rate of 3 x 10~^ Mq yr~^ (for a wind speed of lO'^ km s""'^) 
is in the expected range of a Wolf-Rayet star. The optical afterglow brightens 
unexpectedly around day 2. Its temporal and spectral behaviors before day 2 
are not explained by the wind model. 




Fig. 1. Wind interaction model for the afterglow of GRB 970508 (adopted from Cheva- 
lier & Li [^). Ra dio data are taken from Frail et al. and R-band data from 
Sokolov et al. CI. 



The peak flux F^^ and absorption frequency t'a of GRB 970508 are inferred 
to decrease with time, which is consistent with the wind model but not with 
the standard spherical ISM model in the relativistic regime |Q . Frail et al. |Q 
discussed in depth these problems and suggested jet effects and a transition to 
non-relativistic expansion as possible solutions. They proposed three phases of 
evolution in a constant-density medium: a relativistic jet phase followed by a 
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phase of jet spreading, which is in turn followed by a phase of non-relativistic, 
spherical expansion. The jet spreading reduces the peak flux Fij_^ , and the emis- 
sion expected from the non-relativistic phase is shown to be consistent with the 
radio data, although the complete model has yet to be calculated. The relatively 
slow decline of the optical afterglow for over 100 days was unusual compared 
to other well-studied afterglows. 



4.3 Jets, Winds, or Non-Relativistic Evolution? 

Besides the four sources discussed above, there are half-dozen or so GRBs with 
reasonable multi-frequency coverage of afterglows, particularly at radio frequen- 
cies, which allows for detailed modeling. These include GRBs 980519, 991208, 
991216, 000301C, 000418, 000926, and 010222 as of August 2001. The major- 
ity of these sources show clear steepening in the optical light curves (GRBs 
991216, 000301C, 000926, and 010222), while in others the evidence for Hght 
curve steepening is weaker (GRBs 980519 and 991208) or absent (GRB 000418). 
A variety of models has been proposed for these sources. One common conclu- 
sion is that the simplest, standard model involving a relativistic spherical blast 
wave expanding into a constant-density medium is inadequate, and further re- 
finements are necessary. GRB 000418 was modeled by Berger et al. 0, who 
found that both a spherical wind model and an ISM-interacting jet model fit the 
available radio-to-optical data reasonably well. For the remaining six sources. 



Panaitescu & Kumar |80 81| , p2| constructed semi-analytic jet models for both 
constant-density and wind-type ambient media. They found that ISM-interacting 
jet models can plausibly explain the broadband emission of all sources, and that 
wind-interacting jet models are compatible with the emission of GRBs 991208 
and 991216 and can marginally accommodate the afterglows of GRBs 000301C 
and 010222 but cannot explain the observations of GRBs 980519 and 000926. 
Other authors have modeled these sources individually and, in some cases, come 
to different conclusions. We discuss these sources in turn. 

GRB 980519 is one of three sources that show unusually steep light curve 
decay (with a < —2) but weak or no evidence for breaks (the other two being 
GRBs 980326 and 991208). GhevaHer & Li ||| proposed that the steep decay can 
be explained in a spherical wind model, provided that the electron energy index p 
is close to 3, a value that is higher than normally found in GRB afterglows but is 
within the range found in radio supernovae . Radio data provide some support 
to this interpretation [p3 but do not rule out the jet model, first advocated for 



this source by Sari et al. |105| . Indeed, Jaunsen et al. |^ have recently found 
some evidence for a sharp break in the R-band light curve, which they interpreted 
as due to a jet expanding into a wind-type medium. The relative sparseness 
of the data available around and after the break makes its identification less 
secure. If true, the sharpness of the break would be difficult to understand in 
a jet model, especially if the ambient medium is wind-like |6^. The spherical 
wind model of Chevalier & Li did not take into account of inverse Compton 



scattering, which is important for the parameters they adopted (e.g., |103]) 



14 Li and Chevalier 



Inverse Compton scattering lowers the cooling frequency, and poses a problem 
for the wind model in fitting the X-ray afterglow. 

GRB 991208 is unique in that the radio data are sampled well enough to 
allow for a determination of the evolution of the characteristic frequencies Vg, 
and z/„i and the peak flux i^,/,max at Vrn- The inferred scalings i>g, oc ^-0-i5±o.23^ 
Vm t^i '^io.T g^jjj F^^max OC i-0-47±o.20 Qj compatible with the predictions 
of either the ISM-interacting jet model {v^ oc t~^^^, oc t^^, and i^i/.max oc t^^) 
or the spherical wind model {vg. oc t~^/^, oc t^^/'^, and i^,y,max oc t^^/^) 
but not with those of the spherical ISM model {v^ oc i^m oc t^'^/^, and 
-FV.max oc t"). Galama et al. pointed out a problem with the wind model: 
the rate of optical light curve decline implies p > 3.3 and such a large value 
of p would be inconsistent with the radio-to-optical spectrum at day 7.3. The 
inconsistency motivated Li & Chevalier |Q to seek a model with a non-standard 
energy distribution of electrons. They found that a spherical wind model with a 
broken power-law electron energy distribution can reproduce all data as well as, 
if not better than, the jet model. Interestingly, in the jet model of Panaitescu & 
Kumar | ^l[ | for this source, where the jet dynamics is followed semi-analytically, 
the steep decay in optical is mostly attributed to a steepening in the electron 
energy distribution, as in the wind model for this source, rather than jet effects. 

GRB 000301C was modeled by Li & Chevalier |Q using the same spheri- 
cal wind model with an electron energy distribution that steepens at a certain 
(high) energy as the one applied to GRB 991208. The quality of the overall 
fit to the broad band radio-to-optical data is again comparable to that of an 
ISM-interacting jet model The radio data at 8.48 GHz are particularly well 
sampled for this source. The fact that its temporal decay (with a ~ —1.4 Q) 
is much shallower than that observed in R-band (with an asymptotic value of 
a « —2.7) provides some evidence for a steepening of the electron energy dis- 
tribution at high energies (see also [Q). The steepening in energy distribution 
produces a steepening in the optical light curves, which could mimic the jet ef- 
fects. A difference is that the optical spectrum should steepen with time in the 
former case but not in the latter. Rhoads & Fruchter |^ presented IR-optical- 
UV data at several epochs, showing some evidence for spectral steepening with 
time. These data are fitted reasonably well by the wind model. The spectral 
steepening would be more pronounced in the optical to X-ray regime. Unfortu- 
nately, X-ray observations are not available for this source. Kumar & Panaitescu 
[ |63| attributed the steepening of the R-band light curve to a sudden large drop 
in the density of the ambient medium and Dai & Lu p7| to transition to non- 
relativistic evolution. The predictions of these models on radio emission remain 
to be worked out and compared with observations. 

GRB 991216 is a source well observed in radio, optical and X-rays. The decay 
indexes at these wavelengths (a = -0.82 ± 0.02 at 8.46 GHz, -1.33 ± 0.01 at R- 
band and —1.61 ±0.06 at 2-10 kev) are all different There is some evidence 
for light curve steepening at R-band, although the data at late times arc too 
sparse to tightly constrain the time of transition from one power-law decay to 
another and the decay index after the transition Kl[ . The optical to X-ray data 
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are consistent with an ISM- interacting jet model |Q, but the radio data are not 
Frail et al. attributed the bulk of the radio emission to either the reverse 
shock as in GRB 991023 [104| or a second forward shock. Panaitescu & Kumar 
[ pO| examined the multi-frequency data, and concluded that all of the data can 
be fitted reasonably well using an ISM-interacting jet model, provided that a 
large curvature exists in the electron energy distribution, with p changing from 
~ 1.2 to - 2.1. GRB 991216 is therefore the third source (after GRBs 991208 
and 000301C discussed above; see also [^) for which a steepening in the electron 
energy distribution is proposed, although in all three cases the interpretation is 
model dependent. 

GRB 000926 stands apart from other sources in that its X-ray afterglow 
is observed nearly two weeks after the explosion There is clear evidence 
for steepening in the optical light curves, with the decay index changing from 
a = -1.46 ± 0.11 to -2.38 ± 0.07 around 1.8 ± 0.1 days [|3|j9^,||,|l0§ , pointing 
to a jet explanation. However, the decay after the break is significantly steeper 
than that in X-rays, estimated to be a = — 1.89lQi[g which is difficult to 
explain in the standard jet model. Piro et al. [ pO[ proposed an alternative model 
involving a mildly coUimated blast wave expanding into a dense uniform medium 
(n ~ 3 X 10* cm^'^), with the transition to non-relativistic evolution occurring 
~ 5 days after the explosion. Harrison et al. |Q showed that such a high density 
appears to have difficulties reproducing the radio data. They suggested that the 
standard jet model, either in a constant-density or wind-like medium, can fit the 
broad band data reasonably well, provided that the inverse Compton emission 
contributes significantly to the X-ray emission (see also |129] and 

GRB 010222 is well observed in X-rays, with a decay index a = —1.33 ± 0.04 
that is nearly identical to that in the optical after a break around ^ 0.5 days 
(e.g., |69 ). Before the break, the optical light curve is significantly flatter, with 
a ~ —0.6 to —0.8. The steepening of optical light curves can be interpreted as 
due to jet effects |111,101,22|, provided that the electron energy distribution is 
very fl at (w ith p ~ 1.5) as required by the asymptotic relation p = —a after the 
break [105|. This relation, originally derived for p > 2, may not be applicable 
to the case of p < 2 ||,|||. Masetti et al. fH and in 't Zand et al. §§ favored 
an alternative model, in which the transition to non-relativistic evolution occurs 
rapidly in a very dense medium {n ~ 10^ cm^"^). The inferred p ~ 2.2 is more 
in line with those inferred for other GRBs (see however js^). It is not clear, 
however, whether this high density model can reproduce the early radio detection 
of the source [D . 



5 Discussion and Conclusions 

No single model explains all of the dozen or so relatively well observed GRB 
afterglows. The model that comes closest appears to be the ISM-interacting jet 
model at the present time. It fits particularly well the multi-frequency data of 
GRBs 990123 and 990510. To reproduce the afterglow observations of GRBs 
991208, 991216 and 000301C, a relatively flat electron energy distribution with 
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p < 2 is required; the flat distribution must steepen at high energies to p > 2 
for the total electron energy to remain bound. The model also appears capable 
of describing the data of GRBs 000418 § and 000926, provided that inverse 
Compton scattering contributes significantly to the X-ray emission of the latter 
In addition, the steep decay of the optical light curves of GRBs 980326 and 
980519 has been attributed to jet effects [^05|, and GRB 970508 has been mod- 
eled as a mildly collimated jet making a transition to a spherical non-relativistic 
evolution [^ . 

However, a number of the above sources can be modeled equally well or 
perhaps better by the wind-interacting spherical model. This is particularly true 
for GRB 970508. In the case of GRBs 991208 and 000301C, a non-standard 
electron energy distribution is required as is the case for GRB 991216 (Li & 
Chevalier, in preparation). In addition, the wind model can fit the extensive radio 
data of GRB 980425/SN 1998bw, provided that a late energy injection occurs. 
Furthermore, the optical and X-ray afterglow data of GRB 970228 are compatible 
with the wind model, after the subtraction of plausible supernova emission from 
the optical light curve ||9^,^ . Frail et al. |Q searched for this source in radio for 
the first year, and did not detect any afterglow. The upper limits are of order 100 
/xJy or less at 8.46 GHz where the monitoring is most frequent. This absence of 
radio afterglow is not easy to explain in the standard, ISM-interacting spherical 
model, but is compatible with either a spherical wind model or ISM-interacting 
jet model ||T^, since in both models the peak flux decreases with frequency. The 
jet model may have difficulty with the relatively flat decay of the optical light 
curve (a ~ —1.6) unless the electron energy distribution is very flat and/or the 
asymptotic slope after the jet break has yet to be reached. A relatively complete 
optical/IR data set also exists for GRB 980703 P, pl] , p5[ . The inferred large 
extinction intrinsic to its bright host galaxy, coupled with the sparseness of the 
published X-ray |115| and radio data makes it impossible to constrain models 
firmly. 

The presence or absence of supernova emission in the optical afterglow may 
provide a powerful means for corroborating the ISM or wind model. The ISM- 
interacting sources are not expected to be accompanied by supernovae, and one 
indeed finds no evidence for supernova in the most probable ISM interactors: 
GRBs 990123 and 990510. The wind-interacting sources, on the other hand, 
are expected to be associated with supernovae, and the best wind interactors 
GRBs 980425 and 970508 are probably associated with supernovae. The evi- 
dence of supernova association is strong for GRB 980425, as reviewed by Galama 
and Iwamoto & Nomoto in this book. The evidence is weaker for GRB 970508, 
based mainly on a "shoulder" in the late time light curve at the /c-band |10S|; 
supernova emission is not clearly seen in the R-band ||3^ . In the case of an- 
other possible wind interactor, GRB 970228, evidence for supernova has been 
marshaled by Reichart J93 and Galama et al. the alternative explanation 
involving dust echoes |pOt appears less likely [ p4[ , especially if the GRB pro- 
genitors are Wolf-Rayet stars instead of red supergiants . Association with a 
supernova has also been proposed for other sources, including GRBs 980326 |l§, 
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991208 yj, and possibly 000418 Their afterglow observations appear to 

be compatible with both the jet and wind models. In no source can the presence 
of a supernova be excluded definitively |^,^ . 

The expansion into a wind-type medium and the association with supernova 
would point to a massive star origin for GRBs. The massive star connection is 
indirectly supported by the analysis of optically dark bursts by Reichart | |95[ |, 
which indicates that the majority of GRBs may be tied to giant molecular clouds, 
the sites of massive star formation. The connection may be strengthened by the 
presenc e of p ossible Fe line emission in the X-ray afterglows of GRBs 970508 fs^ , 
970828 m 991216 [!§, and 000214 |], which implies a substantial amount of 
circumburst material. We have discussed GRB 970508 as a wind interactor. 
GRBs 970828 and 000214, are optically dark, which is compatible with the idea 
that they are tied to giant molecular clouds and thus (indirectly) to massive stars. 
Their afterglow data are insufficient to constrain ambient interaction models. 
However, some models for the Fe lines require a large Fe mass and densities 
^10* cm"'^ at r « 10^^ cm [ 123| , 114 |, much larger than the densities inferred in 



afterglow models and than expected in the wind from a massive progenitor star. 
The distribution of the dense gas is asymmetric and it is not present along the 
line of sight to the GRB, which is the region probed by the afterglow emission. 
The strong angular variation of the density is surprising; in the supranova model 
[ |114| , the supernova would be expected to affect the density along the line of 
sight. The required Fe mass would be greatly reduced if the lines come from 
a much denser region closer to the explosion center, as in the collapsar-bubble 
model of Meszaros & Rees Q . 

The absence of supernova emission in the afterglows of GRBs 990123 and 
990510 is consistent with their being ISM-interactors with compact star merger 
progenitors. It therefore appears that there are two types of burst progenitors, 
with GRBs 970508 and 980425 representing the best examples of one type, and 
GRBs 990123 and 990510 the other see also The possible presence of 
supernova emission in the optical afterglows of GRBs 970228, 980326, 991208, 
and 000418, coupled with the fact that their afterglows can be fitted by the 
wind model, plausibly put them in the massive star category. This assignment is 
weakened, however, by the fact that the afterglows of these sources can be fitted 
by the ISM-interacting jet model as well. For the other GRBs with relatively 
well observed afterglows, the progenitor types are even less certain. 

The interpretation of GRB afterglows in terms of ambient interaction models 
is complicated by several major uncertainties. For the jet model, the angular dis- 
tribution of matter and energy inside the jet is not clear and the hydrodynamics 
of the jet-ambient medium interaction remains uncertain. Numerical simulations 
are beginning to address these issues [|5^. An open issue common to both the 
jet and wind models is the shape of the energy distribution of the radiative elec- 
trons. Fitting the afterglow data of several GRBs requires a curvature in the 
energy distribution. How the required curvature comes about is not understood. 
This problem is particularly severe for the cases that demand a flat distribution 
with p < 2. In such cases, the afterglow emission depends sensitively on the way 
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in which the energy distribution steepens at high energies, and the steepening 
has not been treated self-consistently. Another open issue is how the postshock 
energy is distributed among electrons, protons, and magnetic field. In addition, 
the spherical wind model faces an interesting dilemma: the wind model implies a 
massive star progenitor, and a coUimated flow is expected if it must escape from 
the center of a star. To some extent, the spherical model can be justified by the 
slow apparent evolution of a jet in a wind |^,^. It may also be possible that, 
upon passing through the star, the jet becomes uncoUimated due to a sudden 
weakening of the lateral confinement. 

We conclude that there is evidence for two types of GRB afterglows in differ- 
ent environments: a constant-density interstellar medium and a wind of possibly 
Wolf-Rayet star origin. The types are not immediately distinguishable partly 
because, at an age of a few days, the preshock wind density is comparable to 
an interstellar density. At an age of seconds to minutes, the preshock density is 
much higher for the wind case, which could make the types more distinguishable. 
Pair production in the ambient medium could substantially modify the afterglow 
emission. Its effects are expected to be far greater in the wind case than in the 



ISM case [112 7C|,[4|. In addition, the large difference in density for the two cases 
at small distances from the center of explosion affects the dynamics of blast 
wave energization, and thus the prompt emission 1 104 ,1^. While some theoreti- 



cal work has been done on the prompt and early afterglow emission, much more 
is needed. Hopefully, rapid follow-up observations of GRB afterglows, to be en- 
abled by HETE-2 and SWIFT, would put tighter constraints on the ambient 
interaction models and the nature of GRB progenitors. 
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